Generalized synchronization of chaos in autonomous systems 
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We extend the concept of generalized synchronization of chaos, a phenomenon that occurs in 
driven dynamical systems, to the context of autonomous spatiotemporal systems. It means a situa- 
tion where the chaotic state variables in an autonomous system can be synchronized to each other 
but not to a coupling function defined from them. The form of the coupling function is not crucial; 
it may not depend on all the state variables nor it needs to be active for all times for achieving 
generalized synchronization. The procedure is based on the analogy between a response map subject 
to an external drive acting with a probability p and an autonomous system of coupled maps where 
a global interaction between the maps takes place with this same probability. It is shown that, 
under some circumstances, the conditions for stability of generalized synchronized states are equi- 
valent in both types of systems. Our results reveal the existence of similar minimal conditions for 
the emergence of generalized synchronization of chaos in driven and in autonomous spatiotemporal 
systems. 

PACS numbers: 05.45.-a, 05.45.Xt, 05.45.Ra 
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Generalized synchronization of chaos is a common phe- 
nomenon occurring in unidircctionally coupled systems, 
where a distinction can be made between a drive or for- 
cing subsystem and another driven or response subsystem 
[l[ . It arises when a functional relation different from the 
identity is established between the drive and the response 
subsystems This phenomenon has been the subject 
of many theoretical and numerical studie s [3, 3 , [H, |g, @] 
and has been observed experimentally pi. IqL lloTlll|. On 
the other hand, there has been recent interest in the 
investigation of chaotic synchronization and other col- 
lective behaviors emerging in networks of mutually inte- 
racting dynamical elements where no external influences 
are present Il3j |. In particular, the phenomenon of 
complete synchronization, where all the state variables 
converge to a single trajectory in phase space, has been 
widely studied in these autonomous dynamical systems. 
In this paper, we present a procedure that allows to ex- 
tent the concept of generalized synchronization of chaos 
found in driven systems to the context of autonomous 
systems. We mean that, under some circumstances, the 
chaotic state variables in an autonomous dynamical sys- 
tem can be synchronized to each other but not to a cou- 
pling function containing partial information from those 
variables. This phenomenon can be seen as a form of 
collective behavior arising in some specifically designed 
dynamical networks. Our procedure is based on the re- 
ported analogy between a sing le driven map and a system 
of globally coupled maps [14j. This analogy provided an 
explanation of dynamical clustering and of stabil- 

ity of steady states in systems with delayed interactions 



We consider a map driven with a probability p, 



Here we search for minimal conditions for the oc- 
currence of generalized synchronization of chaos in both, 
driven and autonomous systems, by using models of cou- 
pled maps. 



Xt+l 



w(xt,yt), with probability p 
f{x t ), with probability (1 — p), (1) 

yt+i = g(yt), 



where f(xt) and g(yt) describe the dynamics of the driven 
and the drive variables, respectively; and the coupling 
relation between them is chosen to be 



w{x t ,Vt) = (1 - e)f(x t ) + eg(y t ) , 



(2) 



where e is the coupling strength. For the driven chaotic 
dynamics we shall choose in most examples f(xt) = 
6 + ln \xt\, where b is a real parameter and xt S (— oo, oo). 
This logarithmic map exhibits robust chaos, with no pe- 
riodic windows and no separated chaotic bands, on the 
interval b G [-1, 1] 0. 

The linear stability condition for generalized synchro- 
nization is determined by the Lyapunov exponents of 
the two-dimensional system Eq. fl}. These are de- 
fined as A x = liniT— >oo lnL x and A y = limT-»oo ln-k 



v 



where L x and L y are the magnitude of the eigenvalues of 

Y\^ = q 3{xt,yt)] 1 ^ T , and 3(xt,yt) is the Jacobian matrix 
for the system Eq.{T]), calculated along an orbit. A given 
orbit {x tl yt} from t ~ to t = T — 1 can be separated in 
two subsets, according to the source of the Xt variable, 
either coupled or uncoupled, that we respectively denote 
as A = {{x tl y t } : x t = w(x t -i,y t -i)} possessing pT ele- 
ments, and B = {{x t ,y t } ■ x t = f(x t -i)} having (l-p)T 
elements. We get 
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where w-, 



dx 



(1 — e)f'(x), and K is a poly- 
nomial whose terms contain products of w x , e, and 
g'(yt) to be evaluated along time. Then L^ = 



Thus we get 



x t eB 



f'(x t )]V T and i y = [nf=o 1 5'(2/ t )] 1/T . 



p In 1 1 



e + lim — 



in n n 

(4) 

1 T_1 

lim -^ln| ff '(y t )| = A ff , (5) 
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where A g is the Lyapunov exponent of the map g(yt)- 
Generalized synchronization occurs if the Lyapunov ex- 
ponent corresponding to the driven map is negative Q; 
i.e., A x < 0. 

Figure [T] shows the Lyapunov exponents of the driven 
system Eq. (JTJ) with g ^ f ', when the probability p is 
varied. The exponent A y is constant and positive for 
the chosen chaotic drive g. On the other hand, there 
is a definite value of the p at which the exponent A x 
changes its sign, from positive to negative, signaling the 
onset of generalized synchronization and the appearance 
of a contracting direction in the dynamics of the two- 
dimensional map Eq. {!]). 
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FIG. 1: Lyapunov exponents A y and A x for the driven system 
Eq.([T} as a function of p, for g(yt) = 0.5 + In \y t \ and f(x t ) = 
— 0.7 + In \xt\, and e = 0.7. 

Figure [2] shows the orbits of the driven system Eq. {T]) 
for different values of the probability p. Generalized syn- 
chronization in this two-dimensional map system is man- 
ifested by the appearance of a strange attractor for values 
of p above some threshold value. 

When g = /, the condition A x < implies complete 
synchronization, where Xt = yt- In this case we get 



A, 



:pln|l-e| + A/, 



(6) 



where A/ is the Lyapunov exponent of the map /. Fig- 
ure[3]shows the stability boundaries, given by A x = 0, for 
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FIG. 2: Orbits of the driven map Eq.JTJ for g(y t ) = 0.5+ln \y t \ 
and f(xt) = —0.7 + In |a;t|, e = 0.7. (a) p = 0.3 (unsynchro- 
nized). (b) p = 0.9 (generalized synchronization); the insert 
is a magnification of the marked square. 



the completely synchronized states of the system Eq. ([T]) 
on the space of parameters (p, e) for different orbits of 
a drive g(yt)- When g = / complete synchronization 
in a unstable periodic-m orbit of the map /, defined by 
f( m \x n ) =W n and satisfying e x J = ^£ =1 |/'(x n )l > 1, 
where {xi,5;2, . . . ,x m } are the set of consecutive points 
on this orbit, can also be achieved in the system Eq. (JTJ) , 
as shown in Fig. [3] 
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FIG. 3: Boundaries A x = (Eq. @) on the plane (p,e) 
for complete synchronization in the driven map Eq. ([T]) with 
f(x t ) — —0.7 + ln|xt|. Dashed line: g(yt) = {xi = 
-0.855762}; solid line: g(y t ) = -0.7 + ln|?/ t | (this bound- 
ary also corresponds to complete synchronization in the au- 
tonomous system Eqs. (O and ([9])); dotted line: g(yt) = 
{xi = 0.18049, x 2 = -2.41208}. 

On the other hand, if g ^ /, the condition A x < 
corresponds to generalized synchronization, character- 
ized by Xt 7^ yt- Consider, for example, the system 
Eq. ([I]) subject to an intermittently applied, constant 
drive g(y) = C, which reduces to a one-dimensional map 
with a Lyapunov exponent A x depending on the param- 
eters C, e, and p. The region where generalized synchro- 
nization arises on the plane (C, e) for a fixed value of p 
is shown in Fig. [¥] 
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FIG. 4: Region for generalized synchronized states Xt =/= C 
(satisfying < 0) for the system Eq. ([1} subject to a con- 
stant drive g(yt) = C, on the plane (C, e). Horizontal lines 
indicate the values C = 2, C = 0, and C = {xi = -0.855762}. 
The synchronization region for C > also comprises the gen- 
eralized synchronized states ((er) < 10~ 7 ) of the autonomous 
system Eq. © with N = 10 4 , having a coupling function 
given by Eq. (|13|l . In both cases, /(a:) = —0.7 + ln|a;| and 
p = 0.5. 



The dynamics of the single driven map Eq. |T]) can be 
compared with the dynamics of an autonomous system 
of maps that share a global interaction with a probability 

P, Vi, 



(1 - e)f{x\) + eH t (xi : j e Qt), with prob. p, 
f(xl), with prob. (1-p), 

(7) 

where x\ (i = 1,2,..., iV) gives the state of the ith map 
at discrete time t, e is the strength of the coupling to the 
global interaction function H; and Q t is a subset having 
q < N elements of the system that may be chosen at ran- 
dom at each time t. Each map receives the same informa- 
tion from the coupling function H at any t with probabil- 
ity p. When this autonomous system gets synchronized 
at some values of parameters, we have x\ = xt- Thus for 
this same set of parameters, the single driven map subject 
to a forcing that satisfies g(yt) = H{x{ = Xt : j £ Qt) for 
long times should exhibit a synchronized state similar to 
that of the associated autonomous system Eq. ([7]). Thus, 
besides complete synchronization where x\ — x t = H, 
other synchronized states, characterized by x\ = x t ^ H 
and which we call generalized synchronization, should 
also occur in the autonomous system Eq. j7J) for appro- 
priate values of parameters. 

A synchronized state can be characterized by the 
asymptotic time-average (a) (after discarding a number 
of transients) of the instantaneous standard deviations 



<j t of the distribution of map variables x\, defined as 

1/2 
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1 N 
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(8) 



where (xt) is the instantaneous mean of the values x l t , Vi. 
Stable synchronization corresponds to (a) = 0. Here we 
use the numerical criterion (a) < 10~ 7 . 

As an example of complete synchronization, consider 
a partial mean field coupling function 



(9) 



In this case the autonomous system Eq. ([7]) can be 
expressed in vector form as 




-Gt f (xt), with probability p, 

i = s v q / 

with probability (1 —p), 
(10) 

where the iV-dimcnsional vectors x t and f (x t ) have com- 



ponents [x t ]j 



and [f(x t )]j = f(x\), respectively, 



I is the N x N identity matrix, and Gf is an JV x iV 
matrix that at each time t possesses q randomly chosen 
columns that have all their components equal to 1 while 
the remaining N — q columns have all their components 
equal to 0. The case q — N and p = 1 corresponds 
to the usual mean field global coupling [3]. The linear 
stability analysis [l9[ of the complete synchronized state 
f{x\) = f(x t ) = H yields 



(1 - c) + -a* 
q 



< l 



(ii) 



where ctk = <5ofc9 [k = 0, 1 . . . , N — 1) are the set of eigen- 
values of the matrix Gt for any t, with the zero eigenvalue 
having (JV — l)-fold degeneracy. The eigenvector corre- 
sponding to k — is homogeneous. Thus only perturba- 
tions of x 4 along the other eigenvectors may destroy the 
coherence. Thus, condition Eq. (fTTj) with k ^ becomes 



pln|l-e| + A/ <0, 



(12) 



which is the same condition for stability of complete syn- 
chronized states in the driven map, Eq. ([6]), when g = f. 
Thus the boundary that separates the region where com- 
plete synchronization occurs on the space of parameters 
(p, e) for the autonomous system Eq. J7J), with H given 
by Eq. ^ for any value of q, coincides with the sta- 
bility boundary A x = in Fig. [3] for the driven system 
with g = f. However, in contrast with the driven case, 
the unstable periodic orbits of the local map / cannot 
be synchronized in the autonomous system because they 
correspond to unstable synchronized states in this sys- 
tem. 
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FIG. 5: Left vertical axis: Bifurcation diagrams of St(N) 
(black dots) and H (grey dots) vs. e for the autonomous sys- 
tem given by Eqs. © and (H|j) with f(x\) = -0.7 + In \xl\. 
Right vertical axis: (a) vs. e (continuous line). Fixed param- 
eters: p = 0.5, N = 10 4 . 



For other functional forms of the coupling function H 
it is possible to find generalized synchronized states in 
the autonomous system Eq. ([7|). For example, consider 
the coupling function 



H 



(13) 



where q < N , r < N , and define the instantaneous mean 
field of the system as 



S t (N) 



1 N 



(14) 



j'=i 



Figure ((5|) shows the bifurcation diagram of St(N) as 
well as (u) as functions of the coupling parameter e for 
the autonomous system with the coupling function H 
given by Eq (flU)) . The mean field S t (N) is chaotic for all 
values of e. For e > 0.54, we obtain (a) = 0, indicating 
that the system is synchronized in a chaotic state for that 
range of the coupling parameter. Figure © also shows 
the bifurcation diagram of the function H with values of 
q and r chosen such that q/r = 2. Note that H = 2 
on the range where synchronization occurs, otherwise it 
is chaotic. This happens independently of the specific 
values of q and r, as long as q/r = 2. Thus, we have 
a situation where dynamical elements in an autonomous 
system converge to a chaotic synchronized state x\ = Xt , 
while a coupling function that contains partial informa- 
tion about the system reaches a value different from 
that state, i.e. H ^ x%. This is the analogous to the 
phenomenon of generalized synchronization of chaos ob- 
served in a driven system Eq. ([I]) with constant drive 
g = C = 2, as illustrated in Fig. 0] Furthermore, for all 
given values of q and r the autonomous system with H 



described by Eq (|13p exhibits generalized synchronization 
that yields H = q/r on the same region of the plane (C, e) 
in Fig. (j4|) where generalized synchronization is observed 
for the map driven with constant g = C = q/r>0. 

The equivalence between an external drive and a global 
coupling function can be used to predict the emergence 
of generalized synchronization in either system from the 
occurrence of this phenomenon in the other. As an illus- 
tration, consider a parametrically, intermittently driven 
map Eq. {T]) with g = 0. Figure [6] shows the region 
where A x < on the plane (p, e) for this system. The 
analogy external drive-global coupling suggests that an 
autonomous spatiotemporal system Eq. (|7|) having a cou- 
pling function that reached a value H — at synchro- 
nization, should possess the same region of stability for 
this state on the plane (p, e). Consider, for example, the 
coupling function 



H = 



1/2 



(15) 



where St {q) is the partial mean field of q maps randomly 
chosen at each time t. This autonomous system reaches 
generalized synchronization, i.e., a chaotic synchronized 
state with x\ ^ H = for any value of q, on the same 
region of the plane (p, e) as in the driven map Eq. ([1]) 
shown in Fig. [5] 
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FIG. 6: Region of generalized synchronization on the plane 
(p, e) for: (i) the single driven map Eq. (ffj with g — (A x < 
0); and (ii) the autonomous system Eq. Q with N = 10 4 and 
H given by Eq. (|15[) for any value of q ((a) < 10 -7 ). In both 
cases, b = —0.7. 

Note that H does not need to be a function of all the 
maps nor it must be active for all times; what matters for 
synchronization is that all elements in the autonomous 
system share the same minimal information at any time. 
Thus, the nature of the common input (external or en- 
dogenous) being received by the local units is irrelevant 
for synchronization in either driven or autonomous sys- 



5 



terns. In both scenarios, each unit evolves as a driven 
map at the local level. 

Generalized synchronization may also arise in an au- 
tonomous network of coupled oscillators that, in addi- 
tion, share a global coupling function H with probability 
p. As an illustration, consider a one-dimensional lattice, 



(l- e - 7 )/(x<) + 2 



H+l 



(1 



eH(x\ 



[f{x\ +1 ) + j(x\- 1 )} + 
j & Qt), with prob. p; 

7)/(4) + i[/(^ +1 )+/(^r 1 )] ) 

with probability (1 — p); 

(16) 

where 7 is the local coupling parameter. The analogy 
can be established with a lattice of similar coupled maps 
subject to a uniform drive g(yt) with probability p, 



(i- e - 7 )/(xj) + i[/(^ +i ) + /(xr 1 )] 

eg(yt), with probability p; 

(i- 7 )/(xj) + i[/K+ 1 ) + /K- 1 )], 

with probability (1 — p). 



(17) 

Periodic boundary conditions are assumed for both sys- 
tems. Figure [7] shows, for a given example, that the re- 
gion of generalized synchronization is the same for both 
systems on the plane (e,7), when g{yt) = H. 
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FIG. 7: Region for generalized synchronization on the plane 
(e, 7) for both, the autonomous network Eq. (|16[1 with H given 
by Eq. (|13p with q = 200, r = 500; and the driven network 
Eq. (|17|1 with constant g = 0.4. For both systems, f(x) = 
4x(l - x), N = 10 3 , p = 0.5. 

In summary, based on the analogy between a single 
driven map and a globally coupled system of maps, we 
have extended the concept of generalized synchroniza- 
tion of chaos to the context of autonomous dynamical 
systems. It means that there can exist a coupling func- 
tion H containing some information about the elements 
in the autonomous system that reaches a state different 
from the state of those elements when they are chaotically 
synchronized. The functional form of H is not crucial for 
achieving generalized synchronization; what matters is 



the sharing of the same information about the system by 
its elements. By comparing a single map driven with a 
probability p and an autonomous system of maps shar- 
ing a global coupling with this same probability, we have 
shown that the minimal conditions for stability of gener- 
alized synchronized states are equivalent in both types of 
systems. The analogy external drive-global coupling al- 
lows to design a spatiotemporal autonomous system ex- 
hibiting generalized synchronization of chaos from the 
knowledge of the occurrence of this phenomenon in an 
associated single driven map. Extensions of this work 
include the possibility of designing global coupling func- 
tions in autonomous systems to achieve specific behaviors 
or patterns as self-organizing phenomena. 
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